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We propose a simple algebraic method for constructing exact solutions of equations of two- 
dimensional hydrodynamics of an incompressible fluid. The problem reduces to consecutively solving 
three linear partial differential equations for a nonviscous fluid and to solving three linear partial 
differential equations and one first-order ordinary differential equation for a viscous fluid. 

PACS numbers: 



I. INTRODUCTION 



The dynamics of incompressible viscous fluid flows is described by the Navier-Stokes (NS) equation. In the regime of 



large Reynolds numbers (Re), turbulence arises, presenting one of the most important unsolved problems in theoretical 
and mathematical physics. Formally, letting Re tend to infinity (or the viscosity tend to zero), we obtain the Euler 
equations. Whose mathematical analysis is still more difficult than the investigation of the NS equation. And the 
NS equations themselves can be regarded at large values of Re as singularly perturbed Euler equations. The results 
by Kato show that the two-dimensional NS equations are globally defined in C° ([0, oo]; iiP(i? 2 )), and that for 
s > 2 and < T < oo the "weak solution" of the two-dimensional NS equations tends to the solution of the two- 
dimensional Euler equation in C° ([0, T]; H S (R 2 )^ pj. In turn, the three-dimensional NS equation is locally defined 



> 

in . 

in C° ([0,t]; H s (R 2 )), and for s > 5/2 its " weak solutions" are approximated by solutions of the three-dimensional 
Euler equations in C° ([0,r]; £P(i? 3 )) when r is determined by the initial data (by the norm) and by the external 
forces ; ■ The nonviscous limit was intensively studied in [j] , Q . The mathematical investigation of the NS and 
Euler equations is undoubtedly one of the most important problems of mathematical physics. The are now enormously 
many publications devoted to various aspects of these equations. Among the most important results, we can mention 
thee proof of the Hamiltonian character of the two-dimensional Euler equations by Arnold [B] and the investigation 
results on the symplectic structure of these equations Q- In @ and @ the Lax representation was found for the 
two-dimensional Euler equation written in the Euler variables. The Lax representation was later also constructed in 
the Euclidean variables (lfj |. 

In this paper, we describe an astonishingly simple, but effective, method for constructing exact solutions of equations 
of two-dimensional hydrodynamics of an incompressible fluid. This method is applicable to both nonviscous and 
viscous fluids (accordingly, to the Euler and NS equations). The paper has the following structure. In Sec. 2, we 
demonstrate this method for the two-dimensional Euler equations describing nonviscous incompressible fluid flows. 
We give some examples of exact solutions in Sec. 3. In Sec. 4, we generalize the above mentioned method to the case 
of viscous two-dimensional flows. We summarize the presented results in Sec. 5. 
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II. TWO-DIMENSIONAL EULER EQUATIONS 



We consider the flow of a nonviscous incompressible two-dimensional fluid. The two velocity components v x and 
v y are expressed in terms of the stream function ip = ip(t, x, y) 



dip dip 
dy ' dx ' 



(1) 



as a result of which the continuity equation divw = 0. In these variables, the two-dimensional Euler equation assumes 
the form fTl|: 



dAip dip dAiP dtp dAiP 
dt dy dx dx dy 



0, (2) 



where A is the two-dimensional Laplacian. 

Equation ^ is a nonlinear equation that is not one of the so-called integrable equations (despite the existence of 
the Lax representation). Nevertheless, it is possible to develop a procedure for constructing exact solutions for it. 
More precisely,there exist transformations that allow finding exact solutions of Eq. ([2]) that describe nonpotential 
(i.e., eddying) fluid flows proceeding from solutions describing potential flows, the crux of the matter is expressed by 
the following theorem. 

Theorem 1. Let ip(x,y) be a harmonic function in a domain D, i.e., = where A is the two-dimensional 
Laplacian. Let F(t,x,y) be a solution of the overdetcrmincd system of linear differential equations 

AF = kF, 

dF _dipdF dip OF dink p ( 3 ) 
dt dx dy dy dx dt ' 

where k(t) is a function of time. Then the function 

= i> + F, (4) 

satisfies in D Eq. ©, i.e. 

<9AV>i dipx <9Ai/>i dip x dAip t 



dt dy dx dx dy 



0. 



The theorem is proved by direct calculation. In what follows, formulas and ^ are called the "dressing" procedure 
for the D = 2 Euler equations. 

Remark 1. The harmonic function ip is regarded as being dependent on two variables, i.e. ip = ip(x,y). But it can 
also be assumed to depend on t as a parameter. It can be easily verified that the theorem also holds in this case. 
Remark 2. Transformation Q resembles the Darboux transformation (DT) used in the theory of integrable systems 
(l2| . Indeed, the essence of the DT consists in determining a solution \& for the Lax pair with a given "incipient" 
potential (which in turn is a solution of the nonlinear equation under investigation) and subsequently using ^ to find 
new potentials. 

The similarity between the DT and the transformation described above is obvious. Indeed, as an intermediate 
step, we must solve two linear equations © with variable coefficient depending on the harmonic function ip, which 
can be regarded as an "incipient" potential because it satisfies Eq.Q and describes a plane potential flow (this is a 
stationary flow if ip is independent of t). In this case, the new stream function ip\ in Q describes a nonstationary 
eddying flow of fluid. Nevertheless, this is not a Darboux transformation, because nonlinear system ([3]) is not a [L, A] 
pair for Euler equation @. The compatibility condition for Eqs. ([3]) has the form 

dk 

Ad-k9= — F, (5) 
dt ' w 



where 



dip dF dip dF 
dx dy dy dx 
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Using Eqs. Q and the fact that ip is a harmonic fvmction,we can rewrite Eq. ([5]) in a different form, 



d 2 *p 
dx 2 



d 2 i> 
dy 2 



d 2 F fd 2 F d 2 F\ d 2 ip 



dxdy \ dy 2 dx 2 J dxdy 2 dt 



1 dk 



-F. 



(6) 



But system ([3]) should not regarded as the Lax pair for Eqs. (J5)) or ©. The point is that a pair of linear equations 
can be regarded as a Lax pair only if their compatibility condition has the form of a nonlinear equation for the 
potentials. In other words, the "wave function" (whose role in the case under consideration would be played by the 
function F), which is merely an auxiliary expression, must not enter the compatibility condition. The situation we 
deal with is quite different because F is contained in Eq. (|6]) . Of course, we can take a next step in finding the 
compatibility condition for Eqs. (|6]) and ((3]), although no elimination of the function F is possible even then. We do 
not discuss the interesting equation of whether this procedure can be stopped after only finitely many equations are 
written. 

Remark 3. System ([3]) involves an auxiliary time function k = k{t). Its expression is not fixed but is determined by 
the form of the harmonic function ip. In particular, is not compatible for all harmonic functions ip system ([3]) i.e., for 
an arbitrary function ip, there can be no function k(t) such that system ((3|) has a solution. 

We show that this elementary approach unexpectedly turns out to be rather effective for constructing exact solutions 
of the two-dimensional Euler equation. 



III. EXAMPLES OF EXACT SOLUTIONS OF THE EULER EQUATION 

Wi consider a harmonic function of the form ip N — X N A N Y N where 



( I x x 2 



,.N 



( 1 \ 

y 

y 2 
\y N J 



(7) 



A N is a constant {N + I) x (N + 1) matrix with elements a^, and i and k range from zero to N . The expression 
for the matrix A N is found by substituting ip N in the equation Aip N = 0. We present several examples for different 
values of N: 



A 1 



0-00 &01 

aio an 



An 



aio an 
-o 02 



A, 





a 00 


aoi 


a 2 


ao3 




aio 


an 


-3a 30 


ai 3 




-ao2 


-3a 03 








I 


O30 


-ai3 









Aa 



aoo 


a i 


a 2 


ao3 


ao4 \ 


aio 


an 


-3a3o 


ai3 





-a 2 


-3a 03 


-6a 4 








O30 


-a i3 











ao4 











/ 



A 5 = 





aoo 


aoi 


ao 2 


ao3 


ao4 


ao5 


\ 


aio 


an 


-3a 30 


ai3 


5050 


3a 33 
10 




-a 2 


-3a 03 


-6a 4 


-10a 05 












a30 


-ai3 


-10a 50 


033 












ao4 


5a 5 
















I 


O50 


3a 33 
10 














/ 



etc. The number M(N) of free parameters(i.e.,the number of independent elements in the matrices A N ) is determined 
as follows: if N is odd,then M(N) = 2(N + 1), and if TV is even,then M(N) = 2N + 1. Because the stream function 
serves only as an auxiliary means for calculating the velocity field and vanishes under differentiation with respect to 
the spatial variables aoo, we can se t aoo = without loss of generality, which is precisely implied in what follows. The 
other entries should be regarded as functions of time, i.e., a,j = &ij{t). Substituting ip N in system ^ and solving it, 
we extract the function F, after which the corresponding function ipi. N i s calculated using formula Q. 
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We demonstrate this approach for some specific examples. Let N = 1. In this case, an = should be chosen. It is 
convenient to parameterize a01(i) and aio(i) as 



aoi(t) = A{t) sina(t), aio(*) — A(t) cos a(t) , 



where a(t) and A(t) are arbitrary real functions. After a simple calculation, we finally obtain 



n 

ipi = A(t) (ysina(i) + xcosa(t)) + (c, singi + c- cosgi) , (8) 

i=i 



where gi — x cos <fii + ysin^ + fi(t), </>,;, Cj, and are arbitrary constants, and fi(t) is given by the integral formula 



fi(t) = /i(0) + / AA(t) sin - a(t)) . 



In this case, k is simply a negative constant |14j, which is set to minus unity without loss of generality. Expression 
© can be straightforwardly substituted in ((5J to to ensure that we have obtained a solution of the Euler equation 
parameterized by two arbitrary functions A(t) and a(t). we take n = 1, A(t) = /i/cosh 2 i/i, a = const, c[ = 0, c\ = c, 
4>\ = <f), and /i(0) = 0. In this special case, the solution takes the form 



= 2 — (y sm a + x cos a ) + c sm 

cosh i4 



where k = \f—k, k < is introduced. Figure 1-3 shows the velocity fields for the corresponding flow at the instants 
t = 0, t = 1, and i = 50. Here, a = 0, <j) = 7r /4, and the other (dimensional) parameters are set to unity. The resulting 
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FIG. 1: Graph of the velocity field for N = 1, at the instants t = in the domain —10 < x,y < 10. 
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a) tanh i4 



(9) 



solution qualitatively resembles the so-called "exulton", a specific solution of the nonlinear Schrodinger equation [13]. 
An exulton is a rational perturbution that appears on the background of a pereodic wave, increases, and then rapidly 
disappears [l5j . Solution ([9]) behaves similarly, namely, at negative times sufficiently large in absolute value, the flow 
looks qualitatively the same as Fig. 3. As t aspiring to zero at the left, the first term in solution ©begins to dominate, 
which distorts this pattern. The maximum distortion occurs at t = (see Fig. 1). After that, the first term disappears 
© exponentially fast, and we again have original (although somewhat displaced) pattern shown in Fig. 3. 

We now consider the case N — 2. Assuming that aoi, aio, an and ao2 are function of time, we seek F in the form 
F = £(i) sin(a(i)a; + b(t)y + c(t)). Solving Eqs. (j3]) and using function ((U we obtain 



(10) 
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FIG. 2: Graph of the velocity field for N = 1, at the instants t = 1 in the domain —10 < x,y < 10. 
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FIG. 3: Graph of the velocity field for N = 1, at the instants t = 50 in the domain —10 < x, y < 10. 



where £o=const, and the functions a(t), b(t) and c(t) are found from system of ordinary differential equations 

= -a n {t)a{i) - 2a 02 (t)b{t), ^ = a u (<)&(<) - 2a 02 (t)a{t), 

(11) 

^ = -aoi(*)o(t) + a w (t)b(t). 
at 

The functions ciij(t) can now be chosen arbitrarily. At the next step,we substitute them in system and solve it 
for a(t), b(t), and c{t). Finally, substituting the resulting functions in (|T0l) we obtain the desired solution. 
As a simple example, we set an =const and ao2=const. Calculating as described above, we find 

a(t) = ao coshcji, b(t) = — (an coshwi + oj sinh tot) , 



where ao is the integration constant and lu = y/ a n + 4oq 2 . For c(t), it is convenient to set aio(t)b(t) — aoi(t)a{t), 
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which gives c=const and we can choose c = without loss of generality. We note that solution (TUB is not singular, 
and we have kit) — —a 2 {t) — b 2 (t). 

The behavior of the velocity field at t = and t = 50 is shown in Fig. 4, 5, where we have aio = aoi = 0, and 
all the other parameters are set to zero. We note that in contrast to the previos solution, at times t with a large 
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FIG. 4: Graph of the velocity field for N = 2, at the instants t — in the domain —0.5 < x, y < 0.5. 
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FIG. 5: Graph of the velocity field for iV = 2, at the instants i = 50 in the domain —0.5 < x, y < 0.5. 

absolute value, the first term dominates instead of the second. Furthermore, for large values of x and y, the first 
term in expression (ITU1) also dominates. Because the distortion of the structure occurs only in the neighborhood of 
the origin as t — > 0, the small spatial domain —0.5 < x, y < 0.5. is taken in Fig. 2; the indicated changes would not be 
noticeable at the scales chosen in Fig.l. 

Finally, we note that the same method can also be used to construct solution describing potential flows. For 
example, taking N — 3 and k = 0, we can obtain the solution of the Euler equation 

i'i = a iy + a 02 y 2 + a 03 y 3 + a w x + a n xy - 3a 3Q xy 2 + a 13 xy 3 - a Q2 x 2 - 3a 3X 2 y+ 
+a 30 x 3 - a u x 3 y + e ax+b v+ c (Asin(bx — ay — c) + Bcos(bx - ay - c)) , 
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where all coefficients except A and B are arbitrary functions of time. It can be easily shown that Aip\ = 0. Of 
course, two-dimensional potential flows can be described using the powerful method based on introducing a complex 
potential [TTJ] , and using the technique described above therefore seems inconvenient in this case. The advantages of 
our approach for constructing exact solutions in the case of two-dimensional Euler equations become obvious only 
when applied to nonstationary eddying flows of fluid. 



IV. TWO-DIMENSIONAL FLOW OF AN INCOMPRESSIBLE VISCOUS FLUID 



However astonishing, the technique described above can be easily generalized to the case of a two-dimensional flow 
of an incompressible viscous fluid described by the equations 

dt dy dx dx dy 
where v is the kinematic viscosity. Namely, the following theorem holds. 

Theorem 2. Let ip be a harmonic function in a domain D, i.e., Alp = 0. Let F(t,x,y) be a solution of the 
overdetermined system of linear differential equations 



AF = u(t)F, 

dF__dipdF_ dipdF_ 
9< dx dy dy dx 



(13) 



uU(u) = vu 2 , (14) 
>n ipi = ip + . 
- = uA 2 ip v 



where u(t) satisfies the ordinary differential equation 

du(t) 
dt 

and U{u) is an arbitrary function of the argument. Then the function ipi = ip + F in D satisfies Eq. (fl 

dAtpi dipx dAipi djpi dAiP 

dt dy dx dx dy 

As in the case of Theorem 2, Theorem 1 is proved by a straightforward calculation. 

As the simplest example, we consider the "dressing" on the zero background, i.e.,"0 = 0. we choose U(u) — vu. As 
follows from Eq. (|14p . expression u is a constant. Solving system (|13[) . we obtain 

tpx = Ke-^ a2+b2 y (Cx sin(ax + by) + C 2 cos(ax + by)), 

where K , a, b, and C± t 2 are arbitrary constants. Of course, this simple and well-known solution, which we present here 
only to demonstrate the workability of the method. 

We now consider case N = 1. By analogy with expression (|8]) , we can construct a solution of Eq. (fl2|) in the form 

ipi = A(t) (ysma(t) + xcosa(t)) + £i(i) sing + cosg, (15) 

where 

g = R(x cos (p + y sin (p + f{t)) , U(u) = vu, u = —R 2 = const, 

A(t), a(t), and f(t) are arbitrary function, <p is an arbitrary constant, and the functions £i{t), and £a(t) are solutions 
of the system of ordinary differential equations 

= -vR% + R - A(t) sin(0 - a(t))) 6, 

^ = -»R% a® -n(0 - «(*))) a. 

Of course, by analogy with expression (j5]), a superposition of arbitrarily many sine and cosine functions can be 
constructed instead of (fT5j) . This observation is a consequence of the linearity of Eqs. (IT3|) and of the linearity of the 
transformed function ip in the spatial variables. In special case df(t)/dt — A{t) sin(</> — a(t)), we obtain 

tpi = A(t) (ysina(t) +xcosa(t)) + e _!/R2 *(£i(0) sing + ^(0)cosg), 

where f 1,2(0) is the integration constant. Comparing this ([8]) and ((9J, we see that including viscosity, as should be 
expected, leads to an additional exponential factor descibing the dissipation. 
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V. 



CONCLUSION 



The equations of two-dimensional hydrodynamics of an incompressible fluid thus admit a simple algebraic method 
for constructing exact solution. Here, we confined ourself to demonstrating the simplest solution for the cases of 
nonviscous and viscous fluid separately. There is no doubt that the suggested technique can be used to construct 
many much more complicated and physically interesting solutions. 

The most interesting open equation is using the "dressing method" to solve boundary problems. For example, we 
consider Eqs. (|T3l) and ([T4|) . If we solve the two equations Atp = and AF — u(t)F with given identical boundary 
conditions S (it is reasonable to consider the second boundary problems, i.e., the Neumann problem), then the solution 
ipi obtained by the above mentioned method will satisfy the boundary conditions 25. On the other hand, the 
function F must satisfy an additional dynamical equation, which may turn out to be incompatible with solution of 
the given boundary problem. Consequently, the equation should be stated as follows. What is the class of Neumann 
boundary conditions that is compatible with system (fT3")) . (|14l) and the equation Atp = 0? We hope to return to this 
problem in our further publications. 
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